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ABSTRACT 

We use a sample of 121 spectroscopically normal Type la supernovae (SNe la) to show that their 
intrinsic color is correlated with their ejecta velocity, as measured from the blueshift of the Si II 
A6355 feature near maximum brightness, usi n- The SN la sample was originally used by I Wang et al.l 
(2009) to show that the relationship between color excess and peak magnitude, which in the absence 
of intrinsic color differences describes a reddening law, was different for two subsamples split by i>si n 
(defined as "Normal" and "High- Velocity" ) . We verify this result, but find that the two subsamples 
have the same reddening law when extremely reddened events (E(B — V) > 0.35 mag) are excluded. 
We also show that (1) the High- Velocity subsample is offset by —0.06 mag to the red from the Normal 

- U max cumulative distribution functions 



subsample in the (B n 



Vmax)-M y plane, (2) the B 



of the two subsamples have nearly identical shapes, but the High- Velocity subsample is offset by 
~0.07 mag to the red in B max — Knax, and (3) the bluest High- Velocity SNe la are —0.10 mag redder 
than the bluest Normal SNe la. Together, this evidence indicates a difference in intrinsic color for 
the subsamples. Accounting for this intrinsic color difference reduces the scatter in Hubble residuals 
from 0.190 mag to 0.130 mag for SNe la with Ay < 0.7 mag. The scatter can be further reduced to 
0.109 mag by exclusively using SNe la from the Normal subsample. Additionally, this result can at 
least partially explain the anomalously low values of Ry found in large SN la samples. We explain 
the correlation between ejecta velocity and color as increased line blanketing in the High- Velocity 
SNe la, causing them to become redder. We discuss some implications of this result, and stress the 
importance of spectroscopy for future SN la cosmology surveys, with particular focus on the design 
of WFIRST. 

Subject headings: supernovae: general — distance scale — dust, extinction 
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1. INTRODUCTION 

Type la supernovae (SNe la) are very good dis- 
tance indicators after making an empirical c orrection 
based on their li ght- cu rve shape and color (|Phillipsl 
119931 iRiess et alJ Il996f >. The distances are precise 
enough to determin e that the expansion of the Uni- 
verse is accelerating (|Riess et alj|1998t iPerlmutter et al.1 
Il999h . and co nstrain the equat i on-of-state parameter of 



Dark Energy (lAstier et al.l2006tlWoqd-Vasev et"aT] | 2007 
Riess et al.l2007tlHicken et al.ll2009bl:lKessler et al.ll2009 ; 



Amanullah et all I2010D . The current best precision 



is 

0.15 and 0.11 mag in absolute magnitude (after correc- 
tions) when using optical information only and combin- 
ing it with near-infrared data, respectively (Mandel et 
al., in prep.), alt hough most samples have a precision of 
-0.17 mag (e.g., Ifficken et al"1l2009al ). 

Since peak luminosity, light-curve shape, intrin- 
sic color, and most spectral parameters all correlate 
with each other a n d with measured 56 Ni mass (e.g., 
iNugent et all 119951: IStritzinger et al.1 120061 ). measuring 
the light-curve shape of a SN la is adequate to deter- 
mine its peak luminosity. Observing a SN la in multiple 
passbands allows one to estimate its reddening. There- 
fore with only light curves in a few passbands, one can 
precisely measure the distance to a SN la. 
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However, as mentioned above, the precision in mea- 
suring the distance modulus of an individual SN la has 
a floor of 0.15 mag when observing (rest-frame) optical 
passbands. There have been many attempts to further 
reduce the Hubble scatter with a "second parameter" 
that does not strongly correlate with other luminosity in- 
dicators. In particular, since spectra must contain more 
information than light curves, several investigations have 
looked for correlations between spectral parameters and 
both light-curve shape and Hubble residuals, with most 
correlations between spectral parameters and light-curve 
shape not improving the Hubble scatter. However, two 
studies have been able to make a further reduction to 
the Hubble scatte r using the ratios of fluxes at partic- 
ular wavelengths (|Folev et al.l 120081: iBailev et al1l2009h , 
but both require further data to verify the results and a 
theoretical understanding of why particular wavelengths 
h ave such powe r. 

iBenetti et al.l (|2005| ) showed that although a single pa- 
rameter can adequately account for the variation in the 
peak luminosity, light-curve shape, and temperature of 
SNe la, the velocity gradient of the Si II A6355 fea- 
ture, vsi ii, is not well correlated with these other ob- 
servables. The Si II A6355 feature, being relatively 
strong and isolated in maximum-light spectra, is the 
hallmark fe ature of SNe la (f or a review of SN spec- 
troscopy, see iFilippenkoll 1 9971 ) . IBenetti et"aTl (|2005l ) sep- 
arated their sample of SNe la into three groups: high- 
velocity gradient (HVG) objects, consisting of SNe la 
with vsi ii > 70 km s~ x day _1 and Ami^(B) < 1.5 mag, 
low-velocity gradient (LVG) objects, consisting of SNe la 
with vsi ii < 70 km s~ x day ~ 1 and Ami5(_B) < 1.5 mag, 
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and FAINT objects which have Amis (B) > 1.5 mag. 
The HVG and LVG subsamples have similar light-curve 
shape, peak luminos ity, and color demographics. Finally, 
iBenetti et all (|2005| ) found that HVG SNe la tended to 
have a higher velocity for the Si II feature near maximum 
t han the LVG SNe l a. 

IWang et al.l (|2009l) . hereafter I W09l presented an analy- 
sis of 158 spectrosc opically normal SNe la (i.e., "Branch- 
normal" objects; [Branch et al.l |1993|: which excludes 
SNe la similar to SN 1991T; IFilippenko et ail fl992b; 
Jefferv et al.l 119921: SN 1991be: IFilippenko et al.lll992a ; 
Leibundgut et al.l 119931; an d more peculi ar SNe la; e.g. , 
SN 2000cx: ILi et al.l l200ll: SN 2002cx; iLi et al.l 120031: 
SN 2 003fg; iHowell et al.n200l SN 2008ha: iFolev et al.l 
2009), which were separated by their Si II A6355 veloc- 
ity near maximum light, vqi jj. They showed that this 
selection was essentially equivalent to using vgi 11 for the 
SNe la where there was a good spectroscopic sequence. 
SNe la above and below a velocity of ^11,800 kms -1 
were classified as "High- Velocity" and "Normal" SNe la 5 , 
respectively. They looked at the correlation between 
color excess (essentially an offset from maximum-light 
B — V color) and peak absolute magnitude corrected for 
light-curve shape, but not for color. If the color excess 
is truly the result of dust, this relationship is between 
color excess and extinction, and the slope of a line fit to 
the relationship should provide a measurement of the ra- 
tio of total to selective extinction, Ay / ' E{B — V) = Ry. 
The objects within the two subsamples follow a different 
relation between peak absolute magnitude (again, un- 
corrected for color) and color excess, indicating different 
values of Ry for the two subsamples: Ry = 2.36 and 
1.37 for the Normal and High- Velocity subsamples, re- 
spectively. Both of these values, are significantly less 
t han the Milky Way val ue of Ry = 3.1. 

iKasen fc Plewal (|2007f ). hereafter lKP07l presented the- 
oretical light curves and spectra of a SN la model from 
several different viewing angles. The ejecta in the model 
were asymmetric such that different viewing angles had 
different values for vsi 11 and vgi 11 but similar values of 
Amis(-B) and peak magnitude. Since it was the same 
explosion model, all generated light curves and spectra 
come from an event with a single total energy and 56 Ni 
mass. Although the details of the model presented by 
IKP07I may not properly describe SNe la, the generic at- 
tributes of the observables are useful when discussing dif- 
ferences in ejecta velocity for SNe la with similar peak 
luminosity and/or light-curve shape. 

Alt hou gh nei ther study focused on this result, both 
IW09I and lKP0"7l found correlations between f si 11 a nd the 
B — V color of the SN at maximum light. W09 found 
a ^0.1 mag offset between B — V col ors of the Normal 
and High- Velocity subsamples. IKP07I showed that both 
B — V color and vsi 11 depend on the viewing angle of 
their asymmetric two-dimensional model. 

This manuscript examines the correlation of color with 
ejecta velocity. It attempts to explain the physical cause 
of the correlation and examines the implications for 
measuring SN la distances. We suggest that a single 



maximum-light spectrum will significantly improve the 
precision of the measured di stance for a given SN la. In 
Section [31 we re-analyze the IW09I sample, examining in 
detail the relationship between these quantities. We find 
that the intrinsic color of a SN la at maximum bright- 
ness depends on the velocity of the ejecta. We also show 
that a proper treatment of the intrinsic colors of SNe la 
with different ejecta velocities can significantly improve 
the precision of distance measurements. In Section [31 we 
provide a simple theoretical explanation for our results 
as well as examine more complex models. In Section [4l 
we summarize our results and discuss their implications 
for future studies. 

2. EJECTA VELOCITY AND INTRINSIC COLOR 



5 We use "Normal" to describe these objects since W09 used 
this terminology. But we caution that it could be confused with 
"Branch normal," which both Normal and High- Velocity SNe la 
are considered. Calling these objects "Low- Velocity" in the future 
may reduce this potential confusion. 



As discussed above, IW09I presented an analysis of 158 
SNe la that concluded that separating the objects into 
two subsamples based on vsi 11 and treating their extinc- 
tion corrections differently will reduce the scatter of their 
Hubble resid uals. I n this section, we re-analyze the data 
presented by IW09I . showing that their analysis stands. 
However, (1) a different Ry for the two subsamples is 
only necessary when including the reddest events, (2) for 
SNe la with E(B — V) < 0.3 mag, both subsamples have 
reddening laws consistent with Ry = 3.1, (3) there is a 
clear offset in the color of SNe la in the two subsamples, 
(4) separating the objects into two subsamples can par- 
tially account for the low Ry values measured for large 
samples of low-reddening SNe la, and (5) using this dif- 
ference in color, one can significantly reduce the scatter 
of the Hubble resi duals for these SNe la. 

From the entire IW09l sample, we remove SN 2006bt, 
which has spectral and photometric peculiarities that 
result in a peak lumi nosity fainter tha n its li ght-cu rve 
shape would suggest ([Foley et al.ll2010l ). Like IW09l . we 
remove SN 20061f, which has a Milky Way reddening of 
E(B — V) = 0.97 mag and was ~0.7 mag brighter than 
that of a typical SN la after light-curve shape correction. 
To avoid large syst emati c uncertainties for the peak ab- 
solute magnitudes, |W09j excluded all SNe with z < 0.01 
that do not have a Cepheid distances from their final 
analysis. We also remove these objects from our sample. 
After making these cuts, 121 SNe la remain in the final 
sample, with 81 and 40 SNe la being in the Normal and 
High- Velocity subsamples, respectively. 

2.1. Verification of Previous Analysis 

To verify the results of lWOl we exa mine their sample 
of SNe la in the same manner as IW09I . In particular, wo 
use Equation 1 of lW09l : 

Ml„ = M zp +a(Am 15 (B)-lA)+R v E(B-V) host . (1) 

This two-co mponent fit i s similar to the two- component 
fit used by iTrippj ([1998]) and in the SALT ([Guv et al.l 
[20051 l2007h and SiFTO (IConlev et al.ll2008h light-curve 
fi tters. 

IW09l provides all necessary information in their Ta- 
ble 1. Fitting t he da ta, we find the same values for M zp , 
a, and Ry as IW09I : -19.26 mag, 0.77, and 2.36; and 
-19.28 mag, 0.75, and 1.58 for the Normal and High- 
Veloci ty SN e la, respectively. After fixing a = 0.75 (sim- 
ilar to IW09D , we refit the data, finding slightly different 
values for Ry (2.48 ± 0.11 and 1.63 ± 0.08 for Normal 
and High- Velocity SNe la, respectively). These values 
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Fig. 1. — Our re-creation of Figure 4 of W09. The light-curve 
shape corrected peak absolute V brightness as a function of host 
color excess. The blue squares and red circles represent the Normal 
and High- Velocity SNe la, respectively. The best fit lines for each 
subsample correspond to Ry = 2.48 and 1.63 for the Normal and 
High- Velocity SNe la (corresponding to blue-black and red-black 
dashed lines), respectively. The Milky Way reddening law of Ry = 
3.1 is also shown as the dashed line. 



differ very slig htly ( and within the uncertainties) from 
those found b y IW091 Figured] presents our recreation of 
Figure 4 from W09. 

From Figure [T] it is clear that the subsamples have 
different relationships between color excess and peak ab- 
solute magnitude after correcting for light-curve shape. 
Similarly, both subsamples have reddening laws different 
from that of the Milk Way. However, a closer examina- 
tion of the data is required. 

2.2. Excluding the Reddest Supernovae 

As noted by IW09I . the reddening laws for the two 
groups are significantly different even if SNe la with 
E(B — V)host > 0.5 mag are excluded from the fit. How- 
ever, SNe la with E(B — V)i lost — 0.5 mag are highly 
reddened objects (corresponding to 0.8 < Ay < 1.6 mag 
depending on Ry). These SNe la are very rarely seen 
at high redshift and only occasionally observed in the 
local universe. In fact, the only high-z SN la in the 
Constitution set of SNe la (jHicken et all [2009b) with 
E(B — V) > 0.5 m ag is from the Higher- z HST survey 
(jRiess et all l2007t l and has E(B - V) host = 0.52 mag. 
This object is also at z = 0.216, making it the lowest- 
redshift SN la of the .ffST-discovered SNe la. For current 
cosmological analyses, this color limit is not particularly 
useful. 

In Figure^ we again present the data of lW09l except 
now we only show the data with E(B — V)host < 0.35 mag 
(corresponding to Ay < 0.57, 0.87, and 1.09 mag for 
Ry — 1.63, 2.48, and 3.1, respectively). For all further 
analysis, the light-curve shape parameter, a, is fixed at 
0.75. As a check, we have also performed all analysis 
leaving this as a free parameter, and it is always con- 
sistent with 0.75. However, for subsamples with a small 
number of SNe la the uncertainty for the parameter can 
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Fig. 2. — Same as Figure [T] except for SNe la with E(B — 
^Ohost < 0.35 mag. The best-fit lines for the restricted subsamples 
are shown as blue-black and red-black dashed lines and correspond 
to Ry = 2.51 and 2.50 for the Normal and High- Velocity SNe la, 
respectively. The best-fit lines for the full subsamples (as shown 
in Figure [TJ are shown as grey dotted lines and correspond to 
Ry = 2.48 and 1.63 for the Normal and High- Velocity SNe la, 
respectively; however, the line for the Normal subsample is almost 
completely covered by the line fit to the low reddening Normal 
SNe la. A fit to the full sample of SNe la with E(B - V) host < 
0.35 mag is shown as a dot-dashed line and corresponds to Ry = 
2.21. 



be large, and fixing the parameter for all fits makes com- 
paring samples more direct and convenient. For the sub- 
sample with E(B — V r )host < 0.35 mag, the data produces 
best-fit values of (M zp ,R v ) = (-19.26±0.03 mag,2.51± 
0.28) and (-19.41±0.08 mag, 2.50±0.47) for Normal and 
High- Velocity SNe la, respectively. The fit parameters 
for the low-reddening Normal subsample are consistent 
with those founds for the entire Normal subsample. Con- 
versely, the low-reddening High- Velocity subsample has 
values of M zp and Ry that differ at the 1.5 and 1.8c level 
from those found for the entire High- Velocity subsample. 
Although these differences are not statistically signifi- 
cant, the 0.12 mag offset in absolute magnitude and 0.87 
offset in Ry (corresponding to a difference of 0.17 mag in 
Ay for E[B — V) = 0.2 mag) can dramatically affect the 
implied distances of SNe la at cosmological distances. 

We further examine this difference by fitting all data 
with E(B — V)host below a certain ceiling and varying 
that maximum value. The resulting best-fit parameters 
are shown in Figure [3] Naturally, the uncertainties for 
M zp and Ry are much larger for a smaller E{B — V) max . 
This is both because of a smaller number of SNe la in 
each subsample and the smaller lever arm for fitting the 
line. Regardless, the two subsamples always have consis- 
tent values of M zp ; however, they are only consistent at 
the ~1ct level for E(B - V) max < 0.4 mag. (Section EU 
shows that this is likely an offset in color rather than in 
peak absolute magnitude.) The two subsamples also have 
consistent values for Ry for E(B — V) m ?x < .4 mag, but 
deviate significantly for higher values (|W09D . Addition- 
ally, both subsamples are consistent with Ry = 3.1 (at 



4 



Foley & Kasen 



Normal SNe la 
High- Velocity SNe la. 
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Fig. 3. — Best-fit Rv (first panel) and A/ Zp (second panel), as well 
as the number of SNe in each subsample (third panel) as a function 
of maximum E(B — V)host ■ The Normal and High- Velocity SNe la 
are represented by the blue squares and red circles, respectively. 
The dashed line in the first panel represents Ry =3.1. 



lcr) for E(B - V) max < 0.3 mag. 

The fit parameters do not change significantly for any 
subsample of Normal SNe la. But Ry changes signif- 
icantly depending on the maximum reddening for the 
High- Velocity SNe la. It is possible that the highly red- 
dened High- Velocity SNe la have a value of Ry that 
is significantly smaller than that of less reddened High- 
Velocity SNe la or any Normal SN la. A low value for Ry 
has been measured for individual High- Veloc ity SNe la 
(e.g.. lKrisciunas et al.l l2000; Wa ng et al.ll2008h . We note 
that SN 2003cg, which is defined as a Normal object in 
the IW09I sample, but is excluded from our analysis be- 
cause of its redshift, w as measur e d to h ave Rv = 1-80 
(jElias-Rosa et al.ll2006D . iGoobarl (|2008f ) suggested that 
a low value for Ry may be the result of multiple scatter- 
ings off circumstellar dust. In this scenario, it is natural 
that low values of Rv are associated with large values 
of E(B — V). Although the possibility of circumstel- 
lar dust for all High- Velocity SNe la cannot strictly be 
ruled out, the higher value of Rv for the low-reddening 
sample suggests that it is not the dominant component 
of any extinction. Furthermore, as seen in Figure [H 
there is an offset between the Normal and High- Velocity 
subsamples in the My-E(B — T^)host plane that must 
be the result of some linear combination of My and 
E(B — V)\ los t- Therefore, the subsamples have intrin- 



sically different light-curve shape-corrected peak lumi- 
nosities (with the High- Velocity subsample being more 
luminous), have intrinsically different colors (with the 
High- Velocity subsample being redder), and/or intrinsi- 
cally different dust distributions (with the High- Velocity 
subsample having more dust along the line of sight). The 
highly reddened SNe la probably have different proper- 
ties from the low-reddening SNe la, with the highly red- 
dened High- Velocity SNe la having a lower Ry than the 
low-reddening members of that subsample. Further in- 
vestigations are necessary to determine (1) why highly 
reddened High- Velocity SNe la have a low value for Ry, 
and (2) if low-reddening SNe la have a value of Ry that 
is inconsistent with 3.1. 

The results of restricting the samples by reddening 
show that the extremely reddened objects highly influ- 
ence the fit parameters. Depending on choices made to 
cull the sample or (particularly for high-z SNe) selection 
effects, these values may vary dramatically. 

2.3. Color Offsets 

By restricting the lW09l sample to have E{B — V^host < 
0.35 mag, we found in Section |2"T21 that the two velocity- 
selected subsamples have different magnitude zeropoints, 
M zp . However, this is only one way to interpret the data; 
the other possibility is that there is an offset in E(B — 
V r )host for the two subsamples. This is effectively shifting 
the lines in Figure[3]left-right rather than up-down. From 
the best-fit lines, the offset is ~0.06 mag. 

To clearly see this effect, we look again at the data from 
rWOl . Inst ead of using the derived E(B — V)\ ms t values 
from IW09I . we use their values of B max — V max (in most 
cases, there is a simple offset between the two values). 
This pseudo-color is simply the B magnitude at B max mi- 
nus the V magnitude at Vmax, with no correction for host 
reddening. We also restrict the sample to SNe la with 
1 < Ami 5 ( B) < 1 .5 mag and B max - V max < 0.5 mag. 
As noted by IW091 the two subsamples have similar de- 
mographics over this range of /S.m\^(B). 

W09 noted that for B max — V max < 0.2 mag there was 
a 0.08 mag offset (0.10 mag if restricting the /S.m\^{B) 
range) in the average B max — Vm ax between the two sub- 
samples. There is also a lack of very blue (B max — V niax < 
—0.05 mag) SNe la in the High- Velocity subsample. The 
bluest objects in a sample represent the lowest redden- 
ing and intrinsically bluest SNe la for that sample. For 
the Normal subsample to have several SNe la bluer than 
the bluest High- Velocity SN la, all High- Velocity SNe la 
must be significantly reddened and/or the High- Velocity 
and Normal SNe la have intrinsically different colors. 
In the top panel of Figure [5J Figure [T] is remade, ex- 
cept each SN la is plotted as a function of B max — V max 
instead of E(B — V)host and only show SNe la with 
1 < Ami5(-B) < 1.5 mag. In the bottom panel of the 
figure, the same data are presented, except the High- 
Velocity SNe la are shifted by B max — Kn ax = —0.06 mag, 
the value found from the best-fit lines. 

Although by eye the overlap between the Normal and 
High- Velocity objects may not appear to be significantly 
different between the raw and shifted data, it is clear that 
the shifted data are a much better match to the bluest 
and reddest Normal SNe la than the raw data. That is, 
this simple shift adequately accounts for the lack of blue 
objects in the High- Velocity subsample and removes the 
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Fig. 4. — Same as Figure[2] except data are plotted as a function 
of Bmax — Vmax and only SNe la with 1 < Ami^(B) < 1.5 mag arc 
displayed. The top panel shows the raw data, while the bottom 
panel displays the same data with the High- Velocity subsample 
shifted by B max — Vmax = —0.06 mag. 



overall offset in the two populations. 

Another way to examine the effect of shifting the data 
is by examining the cumulative distribution functions 
(CDFs) for the subsamples. The CDFs for both sub- 
samples are presented in Figure [SJ Clearly, there is a 
dearth of blue SNe la in the High- Velocity subsample, 
with 26% of the Normal subsample being bluer than the 
bluest High- Velocity SN la. However, the overall shape 
of the CDFs for the Normal and High- Velocity subsam- 
ples are similar. After shifting the High- Velocity sub- 
sample by —0.07 mag, the CDFs are very similar except 
perhaps for the reddest part of the sample where dust 
extinction must be the dominant source of color differ- 
ences. A simple Kolmogorov-Smirnov (K-S) test shows 
that there is a 0.65% chance that the Normal and High- 
Velocity objects have the same parent population. The 
Normal and shifted High- Velocity subsamples have a K-S 
probability of 99.99%, which does not mean that the sub- 
samples have the same probability (especially since the 
High- Velocity colors were arbitrarily shifted by a fixed 
amount), but does show the outstanding agreement in 
the two subsamples' CDFs after the shift. The value of 
—0.07 mag was chosen to maximize the K-S probability, 
but values of —0.06 and —0.08 mag yield values of 73.07 
and 98.96%, resp ective ly. 

As discussed in IW09l . the offset in maximum-light col- 
ors is likely the result of either different amounts of dust 
or intrinsically different colors. If the offset is primar- 
ily from dust, then the subsamples must have intrin- 
sically different light-curve shape-corrected peak lumi- 
nosities. Additionally, for the CDFs to have the same 
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Fig. 5. — Bmax — Vmax CDF for the Normal (blue solid line), 
raw High- Velocity (red dotted line), and High- Velocity shifted by 
—0.07 mag (red dashed line) subsamples. All SNe la shown here 
have 1 < Ami5(B) < 1.5 mag. 



shape, but be offset in color, the High- Velocity SNe la 
must have an additional amount of dust corresponding 
to E(B — V)host ~ 0.07 mag for all values of B 

max Vniax ■ 

Alternatively, a difference in intrinsic color accounts for 
the offset in the CDFs and does not require a difference 
in the light-curve shape-corrected peak luminosities for 
the two subsamples. The lack of High- Velocity SNe la 
with colors as blue as the bluest Normal SNe la is addi- 
tional evidence that there is almost certainly an intrinsic 
color difference. 

2.4. Properly Measuring Ry for a Sample of 
Supernovae 

Several groups have noted that the Hubble residuals 
of SNe la are minimized by choosing a ve ry low (and 



potentially unphy s ical) value of Ry (e.g., [T ripp 



1998 



2007 



| Astier et ail [20061: iGuv et all l2007t IConlev et al . . 
Kessler et al.l2009t[Hicken et al.ll2009UIAmanullah et all 
20101) . The methods vary in detail, but all are essentially 



equivalent to fitting a line to the data in a plot like Fig- 
ure [1] or [21 If SNe la have intrinsically different colors, 
this method ignores critical information. 

An example of this problem is shown in Figured We 
fit the entire sampl e of SNe la with E(B — V)host < 
0.35 mag from W09 with a single line. Since the Normal 
subsample has bluer objects for the same My than the 
High- Velocity subsample, this process will naturally re- 
duce the slope of the fit — effectively reducing the value 
of Ry . Where both the Normal and High- Velocity sub- 
samples are best-fit with Ry = 2.50 separately, combined 
they are best fit by Ry = 2.21 ± 0.21. 

If one properly separates these groups, the situation 
becomes the same as described in Section 12.21 This 
clearly affects the distances measured in previous cos- 
mological analyses and should be accounted for in the 
future. 

2.5. Improvement to Measured Distances 
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maximum B maK — V max color. The dotted and solid lines represent 
the full sample (with 1 < Amis(B) < 1.5 mag) assuming a single 
intrinsic color and separate intrinsic colors for Normal and High- 
Velocity SNe la, respectively. The blue squares and red circles 
represent the Normal and High- Velocity subsamples, respectively. 
Using separate intrinsic colors for the two subsamples significantly 
reduces the scatter. Additionally, Normal SNe la appear to have 
significantly smaller scatter than High- Velocity SNe la. 
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maximum 

The empty and filled black diamonds represent the full sample 
(with 1 < Amis (5) < 1.5 mag) assuming a single intrinsic color 
(labeled 'Single') and separate intrinsic colors (labeled 'Dual') for 
Normal and High- Velocity objects, respectively. The blue squares 
and red circles represent the Normal (labeled 'N') and High- 
Velocity (labeled 'HV') subsamples, respectively. Only the values 
for -B m ax — V m ax < 0.2, 0.4, 1, and 2 mag (from right to left) for 
each sample and method are plotted. 



IW09I noted that using a different value of Ry for the 
two subsamples would reduce the scatter of Hubble resid- 
uals for their sample. However, as Section 12.21 showed, 
the two subsamples have similar values of Ry if the ex- 
tremely reddened SNe la are excluded, and it is not clear 
that two values of Ry is particularly useful for cosmo- 
logical samples. In this section, we examine how the 
Hubble scatter is affected by separating the SNe la into 
two groups based on their ejecta velocity. 

Figure [6] presents the weighted residual scatter, defined 



E(A/i 2 Ax 2 
E(i/^ 2 



1/2 



(2) 



where A/i and a are the Hubble residual the uncer- 
tainty of the distance modulus for a given SN la, respec- 
tively, for different samples as a function of maximum 
-B max — Vmax color. In the figure, the full sample (with 
1 < Ami5(-B) < 1.5 mag) is represented by the black 
lines, while the Normal and High- Velocity subsamples 
are represented by blue squares and red circles, respec- 
tively. The fits which generate the black dotted line as- 
sume a single intrinsic color for the entire sample, and 
allow M zp and Ry to vary for each maximum value of 
-B max — 7 max . This is the approach that has normally 
been used in cosmological analyses. To determine the 
scatter for the other samples, all parameters were fixed 
for all maximum values of i? ma x — Knaxj an d only a color 
offset between the Normal and High- Velocity SNe la was 
applied. 

For B m ax — Knax < 1 mag, the scatter is consistently 
improved by using two intrinsic colors. For SNe la with 
Bmax — Vmax < 0.2 mag (corresponding to E{B — V) < 
0.3 mag and Ay < 0.7 mag), the scatter decreases from 



0.190 mag to 0.130 mag by adopting this method. Fur- 
thermore, the Normal subsample has a scatter of only 
0.109 mag for this color cut. This is better visualized in 
Figure [JJ where only a subset of the data from Figure [S] 
is plotted. 

Increasing the maximum color increases the scatter 
of the High- Velocity subsample, but the scatter of the 
Normal subsample remains constant for all color ranges. 
This is a consequence of the Normal subsample having a 
single value for Ry for all colors while the High- Velocity 
SNe la have a significantly different value for Ry depend- 
ing on the maximum color (see Section [2~2l and Figure [3]). 
This, in turn, increases the scatter of the full sample for 
large values of i? max — Vmax, even when using different 
intrinsic colors for the subsamples. Perhaps counterintu- 
itively, the scatter of the full sample with a single intrinsic 
color increases as the maximum color is decreases. This 
is mostly because the fraction of Normal to High- Velocity 
SNe la increases with increasing color in our sample. 

Despite the relatively small number of High- Velocity 
SNe in this sample, it is noteworthy that the Normal 
SNe have a significantly better Hubble residual scatter — 
even after treating each subsample separately. A future 
study may find a color parameter that scales with veloc- 
ity, rather than a simple splitting of the sample, will help 
improve this measurement. But in the era of large SN 
samples, where systematic uncertainties are larger than 
statistical uncertainties and the sample sizes are large 
enough to choose the "best" objects, using only Normal 
SNe la for measuring cosmological parameters may be 
a prudent choice. Of course it goes without saying that 
such a selection (at least at this time) would require a 
spectrum. 
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Our analysis is relatively simplistic. Peak absolute 
magnitudes are corrected for light-curve shape using only 
a linear relationship with Ami5(i3), the reddening esti- 
mate is derived essentially only from -B max — V max , and 
the distance moduli are derived from peak absolute mag- 
nitudes in a single band (V). Furthermore, the SN la 
sample is split by velocity rather than having a model 
with veloci ty as a parameter (this is currently necessary 
since IW09I has not released their wg; n measurements) . 
Nonetheless, this initial analysis provides the foundation 
for more sophisticated analyses. In particular, vsm (or 
perhaps a more general ejecta velocity parameter) ap- 
pears to be a true "second parameter" which significantly 
improves the precision with which we can measure dis- 
tances to SNe la beyond that of just light-curve shape 
and color. 

3. THEORETICAL UNDERSTANDING 

Thus far, empirical relationships between ejecta veloc- 
ity and intrinsic color have been presented. In this sec- 
tion, we present a simple explanation for these correla- 
tions. We then suggest that this scenario is a natural pre- 
diction of asymmetric SN explosions. Finally, we show 
that a set of model light curves and spectra also show a 
correlation between color and ejecta velocity. 

3.1. A Toy Model 

An extremely simplistic model of a SN SED near max- 
imum brightness is the following: a blackbody spectrum 
is generated from electron scattering in the ejecta, im- 
printed on this spectrum are broad P-Cygni lines from 
strong lines, and this causes significant absorption at 
wavelengths shorter than —4300 A, where line-blanketing 
from Fe- group elements is significant. The wavelength 
where line opacity and electron-scattering opa city are 
similar is at -4300 A (jKasen fe Wooslevl l2007t ) . which 
is near the peak of the B band. 

As ejecta velocity increases, one expects the width of 
line- forming regions in velocity space to increase (i.e., at 
higher velocities, lines get broader); this is especially true 
for saturated lines. Therefore, individual lines overlap 
more for higher-velocity ejecta. For wavelengths domi- 
nated by electron scattering, this should not significantly 
affect the SED. However, for regions where line opacity is 
dominant, or in regions where the two effects are of simi- 
lar strength, higher velocities should increase the opacity. 
Since the B band is a region where the line and electron- 
scattering opacities are of similar strength, the -B-band 
flux should decrease with higher velocity. However, the 
V band, which is in a region dominated by electron scat- 
tering, will be less affected by higher ejecta velocity, and 
the l/-band flux should not change significantly. There- 
fore, one expects that higher ejecta velocity should result 
in a redder B — V color. 

SN ejecta have a negative temperature gradient with 
radius. Since the outer layers are cooler, SNe with higher 
velocities will have more absorption from lower excitation 
(e.g., Fe II vs. Fe III) lines. Because of the larger number 
of optical lines for the lower excitation states, higher- 
velocity ejecta should produce additional line blanket- 
ing at blue wavelengths. Again, one expects that higher 
ejecta velocity should result in a redder B — V co lor. 

This has been shown in some theoretical models. KP07 
generated light curves and spectra for a single asymmet- 



ric SN la explosion mo del fro m different lines of sight. 
The asymmetry of the IKP07I model is the result of a 
surface detonation, but the final ejecta distribution is 
qualitatively similar to that of general off-center explo- 
sions, resulting either from off -center ignition or d etona- 
tion conditions, or both (e.g.. iMaeda et al1l2010a| ). The 
asymmetry of the model is not (at this point) important 
for the discussion, but it allows an examination of SNe la 
with different ejecta velocities for the same explosion en- 
ergy. There is a clear correlation between viewing angle 
and both wgi n and B — V at maximum light. Although 
this was not explicitly mentioned by KP07, the different 
lines of sight showed that with increasing vsi n the B — V 
color became redder. This is a validation of the simple 
explanation presented above. 

There are several ways to generate different ejecta ve- 
locities for a given luminosity, but this is also a natu- 
ral outcome of an asymmetric explosion, where differ- 
ent viewing angles have similar luminosities but vastly 
different ejecta velocities. Multi-dimensional theoreti- 
cal models of SN la explosions h ave shown tha t they 
can be highly asymmetric (e.g., iReinecke et al.l l2002t 
iGamezo et all l2005t iKuhlen et al.l 120061 ) . Further the- 
oretical models have shown that the asymmetry can sig- 
nificantly affect ob servables depending on viewing angle 
(jKasen et a l. 2009). From this work, it is expected that 
SNe la be asymmetric and that the asymmetry to in- 
crease the dispersion of SN distances, which in turn de- 
creases the precision with which SNe la can constrain 
cosmological parameters. 

Spectropolarimetry of SNe la has shown that most 
SNe la have very low continuum polarization, but can 
have significant line polarization, especiall y in the Si II 
A635 5 fea ture Ce.g.. iLeonard et al.1 l2005t IWang et al.1 
120071 see IWang fe Wheeler! 120081 for a review). Spec- 
tropolarimetry probes asymmetries of the SN in the 
plane of the sky. The data suggest that the global asym- 
metry of a SN la is typically < 10%, but that a large 
plume of Si-rich mate rial with a preferr ed axis may be 
present in all SNe la. IWang et al.l (|2007[ ) also presented 
tantalizing evidence that the light-curve shape of a SN la 
is correlated with the amount of polarization in the Si II 
feature. 

Recently, IMaeda et al.l (|2010bf ) found that there were 
significant offsets from zero velocity for relatively un- 
blended forbidden lines in the nebular spectra of SNe la. 
These lines probe the inner ejecta of the SN, and if the 
explosion were symmetric, one would expect th em to be 
centered at zero velocity. iMaeda et al.l (|2010aft found a 
striking correlation between the velocity offsets of these 
features and the velocity gradient of the Si II A6355 fea- 
ture, vgi ii, at early times. This observation connects 
the outer layers of the ejecta to the explosion near the 
center of the star, and suggests an intrinsically asym- 
metric explos i on m echanism for the majority of SNe la. 
iMaund et al.l (2010) also found a correlation between the 
amount of polarization in the Si II feature and «si n , sug- 
gesting that vsi ii is an excellent probe of asymmetry of 
the outer layers of the SN ejecta. 

Our work also fits within this framework. Intrinsic 
color correlates with vsi n, which correlates with vs\ n- 
Since i>si n does not correlate with luminosity, we require 
two explosions with similar luminosities to have vastly 
different ejecta velocities. That is naturally explained if 
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Fig. 8.— Maximum-light Si II A6355 (circles) and Ca H&K 
(sq uares) velocity as a function of maximum-light B — V color for 
the IKP07I model spectra. Each viewing angle is represented by a 
different color with the mapping shown by the color bar at the top. 
The best-fit lines for each feature are shown as black solid lines. 
The dashed line represents the v gj n offset from the mean of the 
Normal velocity clump that W09 used to differentiate the Normal 
and High- Velocity subsamples. It also roughly separates viewing 
angles into two hemispheres. 

SNe la are asymmetric with all lines of sight having a 
similar luminosity but some lines of sight showing signif- 
icantly different ejecta velocities than others. 

3.2. Model Light Curves & Spectra 

The IKP071 models can be used to examine correla- 
tions between various observables with the hope of un- 
derstanding the physics behind our empirical findings. 
IKP07I already noted that B — V color evolution and «si n 
depend on viewing angle. However, they did not exam- 
ine the maximum-light values of these two quantities and 
their correlation. Figure [5] shows the (B — V) ma , x color 
and v$i ii as a function of viewing angle. The two values 
are highly correlated, and a linear fit produces 

v = (-11.1 - 28.0(5 - TOmax) x 10 3 km s -1 . (3) 

A similar relationship exists for the Ca H&K feature, 
with 

v = (-17.8 - 22.7(5 - V) max ) x 10 3 km s -1 . (4) 

Not surprisingly, the velocity of the Si and Ca features 
are highly correlated in the models. However, unlike the 
Si II A6355 feature, which was fit with a Gaussian to 
determine the minimum of the feature, the Ca H&K fea- 
ture had a very complex profile for some viewing angles, 
and the velocity was simply measured from the mini- 
mum of the profile. Thus the precision of this measure- 
ment is limited by the resolution of the model spectra. 
Real SN la data show complex, non-Gaussian line pro- 
files with multiple absorption minima for the Ca H&K 
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feature (e.g.. iBranch et afll2006l ). This is further com- 
plicated by varying quality for each observed spectrum. 
Future work should determine if Ca H&K correlates with 
color in a manner similar to Si II. 

While our toy model and the models of IKP07I predict 
that -B m ax — ^max monotonically increases with ejecta 
velocity, our current analysis simply splits the sample 
into two groups based on vsm- Future studies should 
determine from the data if a linear relationship is more 
appropriate. This simplification may also be why the 
Hubble scatter is larger for the High- Velocity subsample 
than the Normal subsample. Additionally, the velocity 
used to separate the samples is somewhat arbitrary and 
there could be some blending between the samples, lead- 
ing to an increase in the Hubble residual scatter. A linear 
r elation ship would avoid such a potential bias. 

W09 found that the Normal SNe la had a mean value 
of (ugiii) ~ —10,600 km s -1 , and a value of vsm ~ 
— 11,800 km s" 1 was used to split the sample into the 
two subsamples. Using t he sam e velocity difference from 
the mean to separate the IKP07I model spectra (the lower- 
velocity clump is at vsm ~ — 11, 500 km s , so the sepa- 
ration would be at Usui ~ —12, 700 km s _1 ), the velocity 
separation also cleanly separates the models into a lower- 
velocity clump and a higher-velocity tail. Interestingly, 
this separation also almost perfectly separa tes the view- 
ing angles into two hemispheres. The Kl'()7 model colors 
are also systematically redder than the observations, but 
this trend is typical of other light cu rve models (e.g., 
iKasen et al]|2009t iKromer fc Simil2009f ) and is likely due 
to the neglect of non-local thermodynamic equilibrium 
effects in the radiation transport calculations. 

Though the analysis here was done with a single as- 
pherical model, the same sorts of effects can been seen 
in spherical models when the explosion energy is varied 
for a given 56 Ni mass. We ha ve examined the surv ey of 
parameterized 1-D models of iWooslev et al.l (|2007f ) and 
confirmed that a higher kinetic energy results in a redder 
color at peak. 

4. DISCUSSION & CONCLUSIONS 

We have re- analy zed a large sample of SN la data orig- 
inally presented by IW091 We have verified the results of 
IW09I : after splitting the sample by vsi n, the two sub- 
samples have different relationships between the light- 
curve shape-corrected (but not host-galaxy extinction- 
corrected) peak magnitude and color excess, and that 
difference can be explained with different reddening laws 
for the subsamples. However, we find that this conclusion 
no longer holds when the reddest SNe la are excluded 
from the sample, and in this regime, both subsamples 
can be described by the same reddening law. Further- 
more, the two subsamples are offset from each other in 
the plane of color and absolute magnitude. Since the 
bluest High- Velocity SNe la are significant ly red der than 
the bluest Normal SNe la (also noticed bv lW09t ) and the 
CDFs of the two subsamples have the same shape, but 
are also offset in color, we conclude that the two subsam- 
ples have different intrinsic colors. 

This result has significant implications for SN la cos- 
mology. First, by assuming a single intrinsic color for all 
SNe la, one will naturally find a smaller value for Ry 
than if the subsamples are treated separately. Second, 
assuming a single intrinsic color could significantly bias 
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cosmological results, particularly if the fraction of Nor- 
mal to High- Velocity SNe la changes with redshift. whic h 
may have already been seen in data ()Kessler et aLl feOOQ 1 ). 
Finally, we find that the scatter in the Hubble residuals 
for SNe la with B max — V mllK < 0.2 mag is reduced from 
0.190 mag to 0.130 mag by using two separate intrinsic 
colors for the velocity subsamples. The Hubble scatter 
is further reduced to 0.109 mag by exclusively using the 
Normal subsample. To perform this measurement, a sin- 
gle spectrum within a week of maximum brightness is 
required. 

The correlation between ejecta velocity and color can 
naturally be explained by additional line blanketing for 
the higher-velocity SNe la, affecting the SED in the B 
band, but leaving the V band not as greatly affected. 
We show that a simple asymmetric explosion model can 
reproduce the general trends in the data. 

The relationship between ejecta velocity and intrinsic 
color will be a significant tool for measuring extremely 
precise distances with SNe la. Better theoretical under- 
standing of this relationship and additional observational 
investigations are necessary to fully leverage this finding. 
It is also important to incorporate this information into 
existing SN la distance estimators, which currently de- 
pend solely on light-curve information. We have only in- 
vestigated color near maximum-brightness, and it is im- 
portant to see how ejecta velocity correlates with color at 
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